Towards a physical interpretation for the Stephani Universes. 
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A physicaly reasonable interpretation is provided for the perfect fluid, sphericaly symmetric, 
conformally flat "Stephani Universes". The free parameters of this class of exact solutions are 
determined so that the ideal gas relation p = nksT is identicaly fulflled, while the full equation 
of state of a classical monatomic ideal gas and a matter-radiation mixture holds up to a good 
approximation in a near dust, matter dominated regime. Only the models having spacelike slices 
with positive curvature admit a regular evolution domain that avoids an unphysical singularity. 
In the matter dominated regime these models are dynamicaly and observationaly indistinguishable 
from "standard" FLRW cosmology with a dust source. 



I. INTRODUCTION. 



p. The most general class of non-static, perfect fluid solutions of Einstein's equations that are conformally flat is known 

as the "Stephani Universe" ||]. These mathematicaly simple solutions are characterized (in general) by a shear-free 
but expanding and accelerating 4-velocity, with perfect fluid FLRW spacetimes being the acceleration free subcase. 
The sphericaly symmetric Stephani Universe^ and some of their subcases have been examined in numerous papers 
' (see Q for a comprehensive review). For example, as star models |4|, ||^, and as cosmologies generalizing 

FLRW §, H, §1. Other papers have looked at their global properties and singularities [|lO|,||ll|l, 
the thermodynamics of their fluid source [ p^ , p^ , p8[ |, and, more recently, as models with inhomogcncous pressure 
in which to test recent astronomical data related to supernova observations ||2^. Although these solutions are 
popular because of their mathematical simplicity, they are not considered physicaly meaningful because their fluid 
, source seems quite artificial and excessively restrictive: matter-energy density, p, depends only on t while pressure, p, 
00 ' depends on t and r (in comoving coordinates). Another objection to the Stephani Universes is their incompatibility 
with a barotropic relation [Q, |Q of the form p = p{p). However, the lack of fulfilment of barotropic relations is 
not a strong argument to dismiss a given fluid solution, since the latter might be compatible with more meaningful 
_ non-barotropic equations of state. 
O ■ The thermodynamics of the fluid source of the Stephani Universes is given by the balance equations of energy, 
O^' momentum and entropy per particle, as well as by the particle number conservation and the equilibrium Gibbs 
. equation. Previous papers have shown that the general Stephani Universe, having as source a single component perfect 
fluid, is incompatible with the integrability conditions of the Gibbs equation jl^], so that a " thermodynamical 
scheme" is not possible (or in other words, a mathematicaly consistent and closed set of thermodynamical equations 
of state is incompatible with the fleld equations). However, this restriction does not apply to the sphericaly symmetric 
case, or to any subcase whose isometry group has orbits of dimension 2 or larger, hence mathematicaly simple equations 
of state have been found |p7| satisfying energy conditions, but without any basis on actual physical matter models^. 

Accepting the above mentioned limitations, this paper aims at looking for a reasonable (or at least, less objec- 
tionable) interpretation for the Stephani Universes. Using the Gibbs equation as a deflnition of temperature, T, it 
is posible to determine the four free functions by imposing the constraint: p = nk^^T, where n is particle number 
density and is Boltzmann's constant. However, the fulfilment of the relation between p and p characteristic of 
the equation of state of the ideal gas |^ (together with p = nksT) cannot hold because p = p{t) while p = p{t,r). 
We consider the equations of state of: (A) a single component classical monatomic ideal gas, and (B) a mixture of 
non-relativistic and relativistic ideal gases, where the internal energy of the former has been ignored (hence, p and T 
are the radiation pressure and temperature). By writing up the explicit form of the state variables for the forms of 
the free functions that satisfy p = nksT, we provide explicitly the conditions for a tight approximation to the full 
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equations of state (A) and (B) in terms of the convergence of a series of powers of adimensional products of p and a 
parameter Eq, proportional to the energy ratio kBTQ/mc^, where Tq is a suitable constant temperature value and m 
is the particles mass. This series converges in a matter dominated low temperature regime characterized by eg ^ li 
for models having spacelike slices (hypersurfaces of constant t, orthogonal to the 4-velocity) with constant positive 
curvature (like "closed" FLRW cosmologies). Under this approximation, the evolution of the fluid in the matter 
dominated regime becomes very close to a dust FLRW cosmology, to such a degree that observational parameters 
Hq, fio, 9o become practicaly indistinguishable from those of the FLRW spacetime. The same models for which the 
above mentioned series converges are the only ones admiting an evolution range that avoids unphysical features, such 
as a "Finite Density" singularity and an "Asymptoticaly deSitter boundary" (see to [^), and so allowing for a 
regular initial Cauchy hypersurface. Models whose spacelike slices orthogonal to the 4-velocity have zero or negative 
curvature lack a regular Cauchy hypersurface, while the series governing the approximation to the ideal gas equations 
of state (A) and (B) only converge for comoving observers near the symmetry center r = 0. 

Considering their qualities and defects, the subclass of Stephani Universes derived and discussed in the present 
paper provide a reasonable and theoreticaly consistent interpretation for this simple and well established class of 
exact solutions. 

II. THE SPHERICALY SYMMETRIC STEPHANI UNIVERSES. 

The metric in comoving coordinates for the sphericaly symmetric Stephani Universes is usualy given in such a way 
as to emphasize its resemblance with FLRW cosmologies However, we shall use instead the following alternative 
comoving coordinates: 

ds^ = -U^c^df + L'^ [dr^ + f [dO^ + sm^Odip^)] (la) 

where R = R{t), K = K{t) and K r = K,t/R,t — dK/dR, and the functions F(r), f{r) are given by the three possible 
combinations 

/ = r, F = r/2 

/ = sinr, F = sin(r/2) (2) 
/ = sinhr, _F = sinh(r/2) 

The transformation of the radial coordinate r, relating (l)-(2) with the usual radial coordinate for the Stephani 
Universes, r, is given by r = / df/{l + kr'^/A), where k — 0,±1. The matter tensor associated with (1) is a perfect 
fluid 

T"'' = + (3) 

- (4) 

where /i"'' — c^'^u°'u^ + g""^ and p, p are matter energy density and equilibrium pressure. The time coordinate in (1) 
has been selected so that the expansion scalar Q — u"'^ is given by 

The fluid source of the Stephani Universes is shear-free, hence the remaining nonzero kinematic scalar is the 4- 
acceleration: iia = Ua-i,u^, whose explicit form for (1) is 

■ ^ ^A'' = RK^iif 

U 2 [1 + KF-^] [1 + {K- RK^r) F2] ^ ' 
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A perfect fluid source like (2) satisfies energy-momentum balance: T"'' -i, = 0, particle number conservation {nu'^).a = 
and entropy balance: {nsu°')-a = 0, where n, s are particle number density and entropy per particle. These balance 
and conservation laws take the following forms for (1) and (2) 

WaT"^6 = P+(p + p)e = (7a) 

habT'" ;c = ha 'pfi + (P + P)ua = (76) 

(nu");a = n + ne = (7c) 

(nsu");a = (n + nQ)s + ns = s = {7d) 

where A = u°- A, a = u* A^t is the proper time derivative for any scalar function A and O, iia are given by (5) and (6). 
The state variables must also comply with the equilibrium Gibbs equation 

Tds = d(9+pd(l) (8) 

where T is absolute temperature. 

The field equation G\ = {d,iiG / i(?)T\ for (1), (2) and (3) is 



3c2 ^ \ R I B? 

where: 



f 0, / = r, F = r/2, 
fc = <^ 1, / = sin(r), F = sin(r/2), (10) 
l-l, / = sinh(r), = sinh(r/2), 



while the form of p follows from (7a) as 



Rp,R R 
t = -fj- - - 

The momentum ballance (7b) is 



1 + {K- RK^r) F2 



(11) 



fR^p,RK,n 

Q[l + {K-RK,R)F^f 

where a prime denotes derivative wrt to r and (6) and (11) were used to eliminate iia and p+p. The particle number 
density is given by integrating (7c) 

n = ^[l + KF'f (13) 

where A'' is the conserved particle number distribution. Since the spacial coordinates in (1) are comoving, the entropy 
balance (7d) yields simply s = s{r), and so, with the help of (7a) and (7c), the Gibbs equation (8) reduces to 



Ts'={p + p)(^^^ (14) 



with s' ^ 0, since the fluid of Stephani Universes is not isentropic (or equivalently, does not admit a barotropic relation 
p = p{p)). 
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III. DETERMINATION OF THE FREE PARAMETERS. 



The Stephani Universes contain (in general) two time dependent free functions {R{t), K{t) in (1)) whose determi- 
nation is elusive. The imposition of an equation of state (ie a constraint relating state variables) should lead to the 
determination of these free-parameters, howewer not all relations between the state variables is compatible with these 
free functions and with the radial dependence of the state variables (in particular, the barotropic relation p = p{p) is 
incompatible, see ||2^ ). Conversely, if the two free functions are determined by suitable boundary conditions, or 
prescribed in order to obtain specific kinematic effects, a specific constraint among state variables ensues, though such 
constraint is a "formal" equation of state and might be wholy unphysical. The usual boundary conditions invoked 
for determining the free functions are: a matching with a Schwarzschild vacuum exterior along a comoving boundary 




, or the desire to obtain specific kinematic 
or a metric element that is conformal to that 



and/or a relation between p and p along a symmetry center 
efects, like bouncing [0 or acceleration of comoving observers 

of a FLRW metric [g|7 ffj], |19|, Besides the time dependent free parameters, a proper thermodynamical 

approach produces two extra free functions of r: the initial particle number distribution, N(r), and the entropy per 
particle s(r). 

The Stephani Universes, as exact solutions subjected to the thermodynamical laws (7c), (7d) and (8), become 
determined once the four free functions i?(<), K(t), N{r), s(r) are either prescribed, found or obtained. From equa- 
tions (1)-(14), once p and K are known as functions of R, the field equation (9) becomes the Friedmann equation 
whose integration yields i? as a function of t. However, the functions N{r), s{r) cannot be determined from the 
geometry but from thermodynamical criteria. The usual strategies to determine the free functions have lead, in all 
cases, to equations of state whose physics has still remained unclear. We propose in this paper another strategy for 
the determination of these free functions. 

Consider the following relation associated with the equation of state of an ideal gas: 



p — nkgT (15) 
where is Boltzmann's constant. Using p, n from (11) and (13), the fulfilment of (15) implies the constraint 



which becomes explicitly 

(1 + KF'^) 



, R I 1 



s'N 



RF 



P 
n 



2kBF' 



Kp^R - pK^rF 



^ 



(16a) 



(166) 



The only physicaly reasonable solution of (16) (other solutions yield negative values of n and T) follows if the terms 
in square parenthesis vanish independently. Since p, K depend on R while F^ N, s depend on r, we have 



Rp.B. 
3p 



s'N 



(17a) 



2kBF' ^ K^r/K ^ ^ 
s'F pmIp 

where a, h are arbitrary nonzero real constants. Equations (17) lead to 



3a / \h / 7-> \ 3afc 



(176) 



(18a) 



2(1 — a) akvt 
logTV = ^ ^ ' logF, s = So + — ^ log A^, a ^ 1 



l-a 



So + ^— logF, 





(186) 
(18c) 
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Inserting (18) into (9), (11), (13) and using (15) leads to 



R 



8ttG (Ro 



kic^ ( Ro 



i?2 V 



i?2 



P = 



{a-l)po + [il-3b)a-l]poki{Ro/Rf^''F^ /i?o' 



l + (l + 3a6)fci [Ro/Rf'' F"^ 



R 



(19a) 



(196) 



n 



73 \ 3a6 



i?g (a - l)po + [(1 - 36)a - 1] pofci {Ro/Rf"'' 



N 



1 + fci {Ro/Rf^ F2 1 + (1 + 3a6)fci {Ro/Rf^ F^ 



Ro 



3(a-l) 



(19c) 



(19d) 



Prom (18d), (18d), (19c) and (19d) it is evident that the choice a ^ 1 leads to a diverging temperature along a 

symmetry center marked by F(0) = 0. Hence, wc must consider only a = 1, 6 > 0, but in this case, in order to have 
positive p and T, we must also assume ki < 0. It turns to be convenient to redefine the constant parameters 6 and fci 
as 



6 = 7-1, -fci 



£0 



3(7-1) 

where eq is a positive constant. The metric functions in (1), 4-velocity and 4-acceleration (4) and (6) become 

R 



.r_ ^-£^{l2hl)F\RolRV^ 
l-{eoh,)F^{Ro/Rp^ ' 



l-{eo/7i)F^{Ro/Rp 



^co(R/Ror\f 



where 



71 = 3(7-1), 72 =37-2 



The state variables in (19) take the forms 



P = Po 



Ro 



(20) 



(21) 



(22a) 



n ■ 



I - eo{i2 hi) {p I PoV'^F^ \po 
No f p 



Ri \po 



71 \Po 



k„T. 



poSoRl F^ 



No 



1 - (£o/7i) {p/PoV-' F^ 1 - £0(72/71) {p/PoV-' F^ 



\7-l 



while the Priedmann equation (9) is given by 



7-1 



(226) 



(22c) 



(22d) 



(23) 
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This equation must be integrated in order to obtain R = R{t) and have the models fuly determined. 

So far, we have determined the four free parameters of the Stephani Universe by imposing the constraint (15). From 
(18c) and (22a), the state variables in (22b-d) can be all be expressed as functions of two primary thermodynamical 
variables {p,s), so that rewriting (8) in terms of the latter and using (15), it is straightforward to prove that the 
integrability conditions of this Gibbs equation are fuly satisfied. However, the equations of state of ideal gases 
impose, besides (15), further relations among the remaining state variables |23|. We look at two of these equations of 
state in the following section. 



IV. EQUATIONS OF STATE. 

The fulfilment of the relation p — nksT leads to the consideration of various possible ideal gas equations of state 
for the fluid. However, we have: p = p(t), while the remaining state variables depend on both t and r, hence these 
equations of state will not hold exactly, though the constant parameters 7, Eq, po can be determined in such a way 
that the full equation of state holds under suitable approximations. Two cases will be examined separately below. 



A. The classical monatomic ideal gas 

For this case the equation of state is (15) together with |2|] 

p = m(?n + — n fc^T 



(24) 



As mentioned before, an equation like (24) cannot be satisfied exactly. In order to examine under which conditions it 
can be satisfied approximately, we use (15), (22b) and (22c) to write the right hand side of (24) expanded up to third 
order in Eq 



mc n + —nksT 



uic^Nq f p 



Po 



i?3 



2P0 



Po 



37-2 



2(7-1)^° ' 3(7-l)2i?3 



mc^No 



P 



27-1 



Po, 



(87- 



6(7-1) 



2)2 
—po 



mc^No 



27(7-l)3i?3 



PoJ 



37-2 



(25) 



Since po, ^0 are p, R evaluated at present cosmic time t — to when matter dominated conditions arc assumed to 
prevail, then energy density is overwhelmingly dominated by rest mass energy, and so, looking for the best possible 
approximation of (25) to (24), leads to the following identifications in (25) 



Po = 



1 



Inserting (26) into (25) yields 



mc^n ■ 



-nksT 



4/3 



7 = 



13 

T 



(26) 



(27) 



so that the leading nonzero order term in the series inside the square brackets is quadratic in eg. Comparing (27) 
with (24) shows that (22) and (26) lead to a good approximation to (24) provided the power series inside the square 
brackets in (27) converges. A sufficient condition for this convergence requires that each term (with the exception of 
the zero order term) must satisfy 



-) 

PoJ 



25/3 



F2«(eo)««l, e = 2,3. 



(28) 



The form of the adimensional parameter eg follows by using (22b), (22d) and (26) in order to expand p and /c^T 
around Eq 
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2 ( P_ 

PQ 



2/3 



\ 2/3 /X 4/3 

1 + 3 i^% + 6 ^ + 

Po/ \Po/ 



By looking at the common factor outside the square brackets, it becomes natural to identify 

ksTo 



(296) 



(30) 



where Tq is a suitable constant temperature value. With the help of (22), (26) and (30), the exact form of the state 
variables becomes 



l-hoip/pof'F^ \po 



5/3 



1- leoiRo/ RfF^ [r 



(31a) 



N, 



P 



\Po 



X 2/3 

l-£^l^ F^ 

2 \po 



-\ 3 



No 

i?3 



2 V i? 



(316) 



kefo F^ 



\2/3 



2/3 



while the Friedmann equation (22d) is given by 



3R 



kef a F^ 



l-^{Ro/RfF2 l-l£o{Ro/RfF^ 



2 STrGmNo Eqc^ f Rq 



R 



kc^ 



Rq 
R 

(31c) 
(32) 



The degree in which these state variables provide a good approximation to the equation of state (15) (24) strongly 
depends on the behavior of F, p/po and on the value of Eq in (30), a small number for an enormous range of 
temperatures (if m is a protonic mass, we have Eq ~ 10^"* for as high as Tq « 10^ degrees K.). The fulfilment of (28) 
for specific energy and temperature ranges of the models will be discussed in section VII. 



B. Mixture of matter and radiation 



Consider a mixture of two ideal gases, one non-relativistic (superindex "I" ) and the other ultra-relativistic (su- 
perindex "11"), the equation of state in this case is 



p = mc^n^ + -n^keT^ + STi^fc^r" 



(33a) 



(336) 



If the internal energy of the non-relativistic gas can be neglected: n^ksT^ <C n^^kBT^^ but (at least at a given stage 
of the evolution) the energy density n}^kBT^^ is not negligible in comparison with the rest mass energy mc^ri^, then 
(33) can be approximated by ||2^ 



p = mc^n} + Sn^^fcsT 



(34a) 



p = n^^ksT 



(346) 
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where T = T^^ and the particle number densities independently satisfy conservation laws like (7c), therefore n}, n}^ 
are given by expressions identical to (22c) with A^o replaced by A^q and N^. The approximation (34a-b) can be 
especialy suited for a mixture of photons and barions bacause the ratio of the former to the latter is such a 
large number (« 10*), leading to a reasonable description of cosmological mixtures of non-relativistic matter and 
radiation for energy and temperature ranges prevailing after cosmological nucleosynthesis, including the radiative era 
(recombination and decoupling), and up to the present. 

As in the previous subsection, (34a) cannot hold exactly but can be aproximated expanding its right hand side in 
terms of e^. From the forms of n} and T = T^^from (22c) and (22d), and using (34b) in the form 3n}^kBT = 3p with 
p given by (22b), we obtain 



,2^1 



mc n 



3/90 



Po 



37- 



(7-1)^° ' 3(7-l)2i?3 



27-1 



Po, 



aP^ 



-) 

PoJ 



37-2 



3(7-1)2'^" 27(7-l)3i?3. 
so that, considering matter dominated conditions and aiming at the best possible approximation, we identify 



(35) 



Pa = 



1 4 

=3 ^ 7 = - 

7-1 '3 



(36) 



transforming (35) into 



mc n + -n KbI = p 



2/3 

P \ 77.4 2 



Po/ \Po 



, 4/3 
P \ 7:i6„3 



(37) 



indicating, as in the previous subsection, that a good approximation to (34a) requires the convergence of the power 
series inside the square brackets, leading to a condition very similar to (28) 



?/3 



^ F2«(£o)««1, C = 2,3,.. 
Po/ 



(38) 



The form of Eq follows, as in the previous subsection, by obtaining with the help of (22b), (22d) and (36) the forms 
of p and fc^T expanded in terms of Eq, leading to 



mc^N^SoF'^ ( p 



PO 



4/3 



\ 1/3 / \ 2/3 

1 + 2 (-^j Fho + ^i^-^J F'el + 0{el) 



(39a) 



kRT = 



mc^N^eoF^ f p 



1/3 



X 1/3 

P \ 77.2, 



1 + 5 ^ F^eo + 16 ^ F'ei + 0{ei) 
po) \Po 



2/3 

P \ T^2J1 



(396) 



so that it becomes natural to identify 



Nfksh 



(40) 



where Tq is a suitable constant temperature value for the relativistic component. Using (22), (23), (36) and (40), the 
state variables take the following exact forms 



{NfkB%lRl) F^ ( P\^'^ _ (N^'kBTo/Rl) F2 fR 



l-2eo(p/po)'/'i^2 Vpo 



1 - 2£o (Rq/R) F^ \R 



(41a) 



n = 



Po 



1 - eoF' ^ 
Po 



1/3- 



i?3 



2 / ^0 



(416) 
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kBT- 



1/3 



[1 - eoF2 iRo/R)f [1 - 2eoF2 (i?o/i?)] V R , 



(Rt? 



3i? 



i? 



(41c) 
(42) 



The fulfilment of (38) with eg given by (40) controls how close the state variables aproximate (34a). However in 
the present case the ratio Nq^/N^ must be taken into consideration as well. This will be discussed in section VII. 
Independently of the equations of state presented so far, the evolution range of the models is necessarily restricted by 
singularities and other geometric features. We examine these features in the following section. 

V. SINGULARITIES, REGULARITY DOMAIN, SYMMETRY CENTERS AND FLRW LIMIT 



Looking aX p p and T in equations (22), (31) and (41), pending the determination of i? = R{t) by integrating (32) 
and (42), we can identify the following scalar curvature singularities 



R^O 



(43a) 



37- 



3(7-1) 



--'if 



3(7-1) 



(436) 



The singularity (43a) is marked by a hypersurface t — const, all state variables diverge and is similar to a "big-bang" 
in FLRW spacetimes. The type (43b) has been identified in previous literature as a "finite density" (FD) singularity 
[ p2| , p3| . It is, in general, a hypersurface marked by a curve in the (t,r) or {R,r) plane, though it is an unphysical 
singularity because p, T diverge while p remains finite 







[ 3(7-1) 


.Po. 


FD 


[(37-2)£oi^2 



1/(7-1) 



(44) 



for all comoving observers, except along the center worldline r = where F{0) = and so [p/po]fd diverges. Particle 
number density also remains finite, as can be verified by inserting (44) into (22c). 
Another feature worth remarking is the points in spacetime where 



1 



3(7-1) 



3(7-1) 



= 



(45) 



These points mark a spacetime boundary characterized by an asymptoticaly inflationary or (asymptoticaly de Sitter, 
"ADS") behavior |l^: p + p ^ Q and n — + 0. Although ksT = p/n diverges along (45) (because n vanishes for finite 
p), this is not realy a curvature singularity since all polynomial curvature scalars are bounded. 

Since for the cases we are interested (7 — 4/3, 5/3), the FD singularity is marked by larger values of R/Rq than 
the ADS boundary, and since this singularity is spacelike [|l2| 
the solutions all points of spacetime along the fluid worldlines 



13| | , we will consider as the regular evolution domain of 
ying to the future of the FD singularity. These domain 



can be characterized, for every comoving observer, by either one of the following conditions 



R 

Rq 



> 



(37 - 2)eoF^ 
3(7-1) 



1/(37-3) 



Po 



< 



3(7-1) 



1 1/(7-1) 



(37-2)£oF2 



(46) 



From the material presented in previous sections, it is obvious that the FLRW limit of the Stephani Universes is 
given by £0 = if = and N = TVq, so that iia = p' = s' = 0. Since £0 = imphes p = T = 0, this hmit 
leads to a dust FLRW universe, and so a sort of asymptotic dust FLRW state can be reached as the fluid evolves 
in time. Also, for whatever form of R{t) obtained from integrating (32) and (42), the worldlines marked by r = 
where f{rc) = 0, with / given by (2) and (10), characterize the symmetry centers. Along these worldlines, the metric 
function J'gee = Lf, the radius of the orbits of the rotation group, vanishes regularly. Also, iia, p' vanish, as well 
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the radial gradients of all other geometric or thermodynamical variables should also vanish, hence spacetime appears 
isotropic for observers along symmetry centers. Depending on the functions in (2) and (10), the Stephani Universes 
admit one {r — for fco — 0, —1) or two {r — 0,Tr for fco — 1) symmetry centers. 

The singularities (43a-b) and regular boundary (45) have been examined extensively elsewhere, and so we refer the 
reader to the appropriate references (see ||l^ to jl^l). Considering (46) for the cases k — 0,-1 of F in (2) and (10), 
it is evident that these cases do not allow for the existence of a hypersurface t — const in the regularity domain that 
is also a regular Cauchy surface. In the case k — I we have two symmetry centers marked by r = and r — n, and so 
the full domain of the radial coordinate is < r < tt. Since in this case we have < F = sin(r/2) < 1, it follows from 
(46) that a regular Cauchy hypersurface t = ti can be found in the regularity domain, so that a regular evolution takes 
place in the range U < t < Im, where Im marks the hypersurface of maximal expansion characterized by O = 0. In 
this case, the regular hypersurfaces t — const > ti will be diffeomorphic to the invariant spacelike slices of a "closed" 
FRLW spacetime. 



VI. OBSERVATIONAL FEATURES 



As shown in section IV, the fluid sources of the models approximately satisfy the ideal gas equations of state 
(15), (24) and (34), provided conditions such as (26), (28), (36) and (38) hold. The evolution of the models towards 
the present time (t —* Iq, R Rq) approaches that of a FLRW dust cosmology, characterizing a homogeneous and 
isotropic matter dominated universe. It is not surprising that (15), (24) and (34) approach dust sources, the latter 
being realy a convenient approximation of ideal gases in a low temperature and pressure limits. However, the "nearly 
dust FLRW" approximation of the models must be examined by verifying how the observational parameters approach 
as i — > to those associated with FLRW cosmologies. 

The resemblance of the field equation (9) with the Friedmann equation of a FLRW spacetime, together with the 
fact that p = p{t), Q = 0(0; leads to an invariant definition of observational parameters H = Q/3 = L/L = R^t/R 
and n ~ 8ttG / {3c^ H'^) that is also very similar to their FLRW equivalents. Evaluating along R ~ Rq, we have 

3c^ [71 



Rq 



(47a) 



°~3c2i/2 8ttGpoRI/3c'^ + (eo/ii) - k ^ ' 



71/ (HoRo) 



generalizing the FLRW parameters 



fr2 87TG c2 SnGpo SnGpoRl/Sc^ 

^o-^Po-fc;^, no ^ - s.GpoRl/3c^ - k ^^^^ 

where 71 = 3(7 — 1) and all quantities with tildes will henceforth be FLRW quantities. Since we are interested in 
comparing observations in the Stephani Universes and in FLRW dust cosmologies at i? = i?o, it is useful to expand 
Hq and r^o in terms of Hq, VIq up to first order in Eq 



Hf) w Hn 



liiHoRo)' 



;£o + Oiel) 



fin W iln 



MHoRo)' 



r^o + 0{el) 



(HoRoy 



(49) 



showing that, in general Hq > Hq and flo < ^q, and that Hq, fig a-re roughly equal to Hq, Qq plus corrections of 
the order « eo^^o ~ £o- Another important observational quantity is the deceleration parameter q, defined for FLRW 
cosmologies as 



go 



R,ttR 



2 " 
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We can compute, with the help of H given by (9) for the Stephani Universes, an equivalent quantity qo and expand 
it in terms of qq 



Qo 



H.t + 



1 SnGpoRl/Sc^ 



no 2 + 60 



So 



2 87rGpoi?o/3c^ + £o/7i - fc 2 2 - Eq 2 - Eq 



(50a) 



qo~qo + 



no 

71 



{HoRoY 



■so + 0{el) 



(506) 



However the definition of q, unlike H, fl, is coordinate dependent. In order to derive the appropriate equivalent of qo 
for the models based on the Stephani Universes, we need to derive the general form of the relation between red shift 
and luminosity distance for local observations along t = to, r — ro- 

Consider the simplest form of a radial null vector = (fc*, k"^) for the metric (1) with U and L given by (20) 



cdt L 



R 



dv 



U l-eo72J^2(i?Q/i?)3(7-i)^ 



dv 



rk'' = -L 



(51) 



The effects of inhomogeneity and anisotropy on red shift observations 'm t — to of nearby sources (z ^ 1) can be 
appreciated from the relation between luminosity distance vs red shift given for the radial light rays (51) by [p7| 



Ua;bk k 



A4 + ^ 
2 



(52a) 



k'\ k* = -U'^ = 



¥ = tL, Mo = - [ubk^] p AA = AA 



(526) 



where AA is the geodesic parameter of an observed gallaxy. For an expanding, accelerating but shear-free 4-velocity, 
the quantities involved in (52) are 



'^a;bk k — — Uak 



(53a) 



B - 2B^ - • ~ - - 

— fc'^ + Qiiak" - Ufc;cA:^r + UaUbk''k'' 
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(536) 



Evaluating these expressions for (20), (21), (48b) and (47), expanding up to first order in Eq and evaluating at {to, ro) 
yields 



Ho 



1 T -^^eo 
2HoRo 



H 

A£o + ^ <i 3 + 90 + eo 



'^n + „ , (7i+5)/oc (1^2fcfoV 
+ + 2HoRo ^ 2{HoRo? 



Ml + 0{Ml) (54) 



where Ho, qo are given by (47a) and (50). As shown by (54), the red shift observed from a given source depends 
in general on the position of the observer (its radial coordinate ro determining fo,Fo). This is expected, since the 
metric given by (1) and (20) has an explicit radial dependence contained in the functions / = AFF', and so the red 
shift distribution is anisotropic because of the inhomogeneity of the cosmological model. In particular, an accelerating 
and shear- free fluid leads to a monopole plus dipole term in (53a), since the direction of the 4-acceleration (d/dr) is 
a privileged direction. Hence, for a connecting vector in the radial direction A£o has a maximal value, while along 
other directions it must be multiplied by a factor cos^, where ^ is the telescopic angle [^. Observations along the 
symmetry center r — detect an isotropic distribution 



All + OiAel) (55) 

Equation (54) provides then the maximal magnitude of the corrections to local red shift isotropy that inhomogeneity 
of spacetime introduces into to the red shift vs luminosity distance. From (47c) and since /o ~ -fo, these corrections 
are of the order £or2o-Fo ~ £o^^o-Fo ~ ^o^o- For these corrections to be small, we need £o^o ^ Ij which for the models 
considered in sections HI and IV, means that either k — 1 (where < r < n and /, F are finite), or if = 0, —1, the 



Ho 



Z\r=0 = —Mo + -4 

c 2c^ 



3 + go T £0 



2{HoRof 
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domain r of f,F must be necessarily restricted. However, as discussed in section V, the cases fc = 0, — 1 do not allow 
for regular hypersurfaces t = const nor a regular initial Cauchy surface at t = to. For the case fc = 1, if eo 1 the 
corrections due to inhomogeneity are small for all tq compared with the isotropic zero order terms Hq, qq, which means 
that local observations are overwhelmingly determined by Hq, qq, leading to a distribution of red shifts that is almost 
isotropic but whose observational parameters slightly differ from the FLRW parameters Hq, qq. However, because of 
(49) and (50), Hq, qq differ from Hq, qq by corrections of order Eq, and so for sufficiently small Eq, compatible with 
(28) and (38), this difference should be much smaller than the 10-15 percent or so observational uncertainty |2|] in 
the most accurate possible local measurements of Hq- This will be discussed in the following section. 



VII. DISCUSSION AND CONCLUSSION. 



We have presented in previous sections important general features of the subclass of Stephani Universes derived in 
sections HI and IV. We still need to prove that the approximation to ideal gas equations of state (conditions (28) and 
(38)), as well as the approximation to isotropic local observations, are actualy met for physicaly reasonable numerical 
values of the parameters involved. We also need to integrate the Friedmann equations (32) and (42) in order to obtain 
a fuly determined dynamical evolution of the models, allowing us to verify their prediction of the "age of the universe 
today" in comparison with the prediction of standard cosmology (dust FLRW universes). Each case, 7 = 4/3 and 
7 = 5/3, is examined separately below. 



A. The classical monatomic ideal gas source. 



For this model the state variables p,p,n,T are given by (22a) and (31a-c), satisfying (15) and providing (if (28) 
holds) a close approximation to (24). Consider matter dominated conditions along t = to, so that po is matter-energy 
density today. In order to examine the fulfilment of (28), suppose that m is a proton mass and (say) Tq w 100 degrees 
K, leading to Eq w 10^^" in (30). The terms of the form (28) in the series in (27) have the form 



controling the goodness of the approximation of (27) to (24) along the directions of increasing r (increasing F) and 
t (decreasing p). For the cases k — 0,-1, the function F increses monotonously as r increases, so that even for very 
small values of Eq and p/po we might have a diverging series in (27), since F^ < for r > 2. Hence in the cases 

k = 0, —1, a necessary (but not sufficient) condition for (27) to converge into (24) is r < 2, restricting the domain of 
regularity to comoving observers near the center worldline r — 0. In the case k = 1, the function F is bounded and 
satisfies F^ > for all values of (, and < r < tt. Therefore, in this case the convergence of the series in (27) 

and the fulfilment of (56) only depends on the comparison of eq with ratio p/ po for the leading term of nonzero order 
(^ = 2) in (27). From (46) and using the value eo « 10"-^°, the allowed values of p/po (with the maximal value at the 
FD density) is given by 



1 3/2 



15 



< 



PO 



3eoF^ 



10 

sin3(r/2) 



(57a) 



Considering (56) with ^ = 2 for k ~ 1, the conditions for a good approximation of (24) become 



Elsin\r/2) 



PO 



4/3 



10 



-20 



4/3 



< 1, 



10" 



15 P_ 
PO 



< 1 



(576) 



Thus, comparing (57a) and (57b), we see that the approximation to (24) is a good one as long as one excludes very 
early evolution stages close to the FD density. However, from (29b) and (31c), a ratio p/po ~ 10^^ can be associated 
with to T/Tq w 10^", corresponding to relativistic temperatures far above the range validity of the classical equation 
of state (15), (24). Hence, the condition (57b) is perfectly reasonable and allows for a large and physicaly justified 
range of validity for the approximation of (27) to (24), based on the state variables (31). A convenient limit of this 
validity range can be set for p/po < 10^, corresponding approximately to T/Tq < 10**, roughly (and shghtly larger 
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than) the range of post-decouphng temperatures for which the evolution a classical ideal gas without interacting with 
radiation can be justified. The reasoning apphed above for the case fc = 1 is vahd for A: = 0, —1, as fong as we consider 
only spacetime sections with r < 2. 

The range of validity of the ideal gas equation of state affects also the question, raised in the previous section, on 
whether local observations on the models are significantly different from those performed in dust FLRW cosmologies. 
As far as the parameters Ho, fto, qo are concerned, the small value Eq « 10^^" is sufficient to make them (via (49), 
(50b)) practicaly identical to the equivalent FLRW parameters. However, the red shift vs luminosity distance (54) 
contains products of the form eo-f^o ^^"^ ^o/o : and so for the cases k = 0,-1, these terms can be large even for such 
a small Eq. Hence, the same reasoning applies: for observations to be close (up to order eo) to FLRW spacetimes, 
either we assume k = 1 or, if /c = 0, — 1, then we must only consider comoving observers near r = 0. 

The dynamical evolution of the models, especificaly the model k — 1, follows from the integral of the Friedmann 
equation (32), leading to 



Ho{t-U) 



AdA 



[eo/2 + ^ioA - A2] 



1/2 



/iQ arcsm 



[eo/2 + /xoA - A^] 



211/2 



(58a) 



where 



2(f^o - 1) 2(2go - 1) 



£o 



Mo 



SttG 
"3^ 



-PoRq 



(2"£o)f^o (2-eo)'7o + £o 



2(r!o-i) 



2go- 1 



(586) 



and Ai is the value oi A associated with a suitable initial surface marked hy t — ti, for example the surface characterized 
by T/Tq = 10^ corresponding (from (31c)) approximately to R/Ro ~ 10^^. The regularity domain, starting (say) at 
Ai = Ri/Ro = 10~2, must exclude unphysical features like the FD singularity and ADS boundary, marked by the 
coordinate values 



R_ 

Rq 



FD singularity 



(59a) 



_R 

Ro 



= ^Aos = Jf sin(r/2)< J|, 



ADS boundary 



(595) 



corresponding to very early stages of the evolution {R/Rq w 10^'°^), well outside the stage where the initial value 
surface R/Rq — 10~^ has been suggested. This surface is a regular Cauchy surface and the evolution range of the 
models is given by 



A,^ « A, = 10-2 <A< A^,,^ = ^ + « MO (60) 

where ^^^x niarks the surface of "maximal expansion" characterized hy H — 0/3 ~ 0, after which the fluid recoUapses 
following a time symmetric pattern. From (58b) and (60), and bearing in mind that Eq ~ lO^-'^^, we have ^„ax ~ 
ilo/{flo — 1), a very close fit to the value of A^^^ obtained for a FLRW dust universe with fc = 1. 

The time required for the fluid to expand up to A = 1 (the "age of the universe today") and up to A^^^ follows 
directly from (58) and can be expressed as a function of either {Qq,so) or (go,£o)- As shown by figure 1, the curve 
relating i?o(^o — U) obtained from (58a-b) vs fio > 1 for Eq = IQ-^" is almost identical with the curve relating these 
same quantities for a FLRW dust universe with fc = 1. As shown in Figure 2, for any given ilo > 1 the value of 
Ho{to — ti) obtained from a dust FLRW universe with fc = 1 is slightly larger than this same quantity obtained from 
(58). In fact, the ratio of the difference of these values of Ho{to — ti) to the FLRW value is of the order of 10^** of the 
latter value. Therefore, the age of the universe today predicted by this model is practicaly the same as that predicted 
by a dust FLRW universe with fc = 1, with -ffo(^o — ii) ^ 2/3 as ^ 1, thus sharing the same "age problem" 
associated with standard cosmology. Depending on the value of Hq, in the range 50 — 70km(mpc s)"^, standard 
cosmology leads to a range of age values to — ti ^ 10^"^^ years that barely meets accepted age values obtained from 
globular clusters and from nuclear and geological estimates The same result follows for models with fc = — 1 

(assuming r ^ 1), the "age of the universe" is almost the same as that predicted by standard cosmology with fc = — 1. 
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The case fc = (also assuming r <C 1) is slightly different, since in the FLRW models k = necessarily implies 
fig = 1, but for the k = Stephani case the value fc = corresponds to fio < 1- We can only have JIq = 1 in the 
Stephani models if fc = 1 and Eq = 3(7 — 1), a possible but undesirable combination of values since £0 ~ 10"^" while 
7 = 4/3, 5/3. 

From equations (31) for fc = 1 we have p{t,0) ~ T{t,0) = along the symmetry center r — 0, while both p and T 
have their maximal value along any surface t — const along the other symmetry center r — tt (see subsection C for 
discussion on these points), that is, there is a "cold center" along r — and a "hot center" along r — tt. However, if 
(28) and (56) hold and for matter dominated conditions, the variation of p and T at any surface t = const is minimal 
(of the order Eq w 10~^") as illustrated by the approximated forms of state variables at i = to 



Po = ^^^, Po^^^^sm\r/2)[1 + O{eo)], no « ^ [1 + 0(£o)] , To ^ fosm\r/2) [1 + O{eo)] (61) 

clearly indicating a matter dominated regime po/po ~ £0 ^ li with very low pressure and temperature values 
providing an excelent approximation to dust-like conditions. The models with 7 = 5/3 might not be realistic becouse 
they do not allow for the description of another gas (CBMW radiation or other rclique gases) as parts of a decoupled 
mixture accompanying the monatomic ideal gas. However, if the matter content of the present universe is dominated 
overwhelmingly by non-barionic CDM, and if we can justify this CDM to satisfy the equation of state of a classical 
monatomic ideal gas ((15) and (24)), then the case 7 = 5/3 could provide a nice toy model of a CDM universe in its 
latter stages approaching near dust conditions. Obviuosly, it would be necessary in this case to infere the mass of the 
particles forming the CDM ideal gas, so that an appropriate value of Eq be given. 



B. The matter-radiation mixture. 



The state variables for the case 7 = 4/3 follow from (41a-c). This case is perhaps more interesting than the 
case examined before because of the possibility of accomodating a joint description of non-relativistic matter and a 
relativistic relic gas. Considering the non-relativistic and ultra-relativistic particles to be barions and photons, the 
ratio Nq^/N^ can be associated with the ratio of photons to baryons, and so we have N^^/Nq w 10^ (considering 
/i « .6 and fls ~ ^0 |29|] ). Assuming m to be a proton mass and Tq ~ 3 degrees K., a good numerical estimate in 
(40) is £0 ~ 3 X 10^^, a value five orders of magnitude larger than £0 in (30). However £0 is still sufficiently small 
so that (38) might hold and (37) might provide a good approximation for (34) over a wide temperature and density 
ranges. Since T under (34a-b) is the temperature of the photon gas (so that To « 3 degrees K), we can set again po 
to be matter-energy density today (overwhelmingly rest mass density) . 

The arguments discussed in the previous subsection regarding the convergence of the series in (27) for cases fc = 0, — 1 
apply also for the series in (37). Hence, we consider only the case fc = 1, though the results can be applied to the 
cases fc = 0, — 1 for comoving observers near r = 0. As before, we use look at (46), now with the value £0 « 10~^, 
obtaining the allowed values of p/ po (with the maximal value at the FD density) given by 
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10 

sinS(r/2) 



Considering (38) with ^ = 2 for fc = 1, the conditions for a good approximation to (34) become 
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Po 



(62a) 



(626) 



Hence, from (62a) and (62b), a good approximation of (37) to (34a) requires the exclusion of evolution stages close 
to the FD density. From (39b) and (41c), the ratio p/po ~ 10^^ corresponds to T/Tq w 10^, a higher value than 
temperatures of the matter-radiation decoupling era, thus the range of validity for the approximation of (37) to (34) 
exceeds that of a decoupled mixture. 

The dynamical evolution for the case fc = 1 follows from the following integral of the Fricdmann equation (42) 
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where: 



6 = 



The regularity domam, startmg at p/po ~ 10^^ or R/Rq « 10 ^ 
boundary marked by 



flo — Eq 
2(r!o-l) 

a £o rnust exclude the FD singularity and ADS 



R 



Rn 



Aj,^ = 2£osin^(r/2) < 2eo, FD singularity 



(64a) 



-B~ = ^ADs = £osin^(r/2) < eo, ADS boundary (645) 

Rn 

L " -I ADS 

indicating that we can set Ai = A{ti) to a value greater than 2eo ~ 10~^ that would correspond to the decoupling 
era [T/Tq « 10^), for example: Ai — . The regular evolution range of the model extends from Ai up to the value 
^j^^^^, marking the maximal expansion = 0, where the fluid bounces and begins to collapse in a time symmetric 
pattern. Evaluating A^^^ from (63) yields 



^MAx = = 7^ ~ (65) 

a value almost equal to that obtained for a dust FLRW universe with fc = 1. The "age of the universe today" can 
be evaluated by setting A — 1 in (63). Ploting the resulting value H(){t() — ti) vs i7o > 1 and comparing with the 
FLRW spacetime leads to curves almost identical to those shown in figure 1. However, the ratio of differences of 
Ho{to — ti) to the FLRW value yields an even closer fit to the FLRW expression (see figure 2). Hence, as in the case 
7 = 5/3, this model presents the same "age problem" of standard cosmology mentioned before. If we assume that 
the non-relativistic component includes cold dark matter (assuming also that it can be described as an ideal gas), we 
would have a value Nq larger by two orders of magnitude, modifying (40) to £o ~ 10^^. However, this alternative 
value of eo does not change the close fit to the dust FLRW universe, in fact, it makes this fit even closer, since the 
only way to obtain values of i?o(^o — U) significantly different from the dust FLRW values is to have a larger Eq « 1, 
though the resulting values of i?o(^o ~ U) are then significantly smaller than those of the FLRW model. As in the 
case 7 = 5/3, this hopeless situation also occurs for the cases with k = 0, —1. 

As in the previous case, from (39) and (41), we also have a "cold" symmetry center r = where p{t, 0) — T{t, 0) = 
(see subsection C for discussion on this aspect), and a "hot" center at r = tt where both p and T have their maximal 
value along any surface t = const. At t = to the state variables take the approximated forms 



Po = ^^, po^^^^^sin'ir/2)[l + 0{so)], no ^ + ^"^ [l + Q(£o)], To « To sin2(r/2) [1 + O(eo)] 

(66) 

so that matter dominated dust-like conditions do prevail ai t — to, with the variation of p and T of the order 
So ~ 10~^. Also, as with the case 7 = 5/3, as long as Fo and /o in (54) are bonded, observation parameters are almost 
identical to their FLRW equivalents. However, the pressure and temperature in this case correspond to the pressure 
and temperature of the MWCB (with 2o = 3 degrees K.). Hence, the ratio Po/po ~ £0 ~ 10^^ from (66), is roughly 
the accepted value of ratio of the energy density of the CMBR (« 3po) to that of visible non-relativistic matter. 
Regarding the temperature, a deviation of order eg « 10~^ gives an excellent fit with the observed deviations of the 
isotropy of the temperature of the MWCB [E^ , |2^ . We might consider including CDM in the matter content of the 
non-relativistic component so that eo ~ 10" .This might sound reasonable, but CDM is most possibly non-barionic, 
and so, as in the case 7 = 5/3, we would have to justify that it satisfies an equation of state like (34). With or without 
CDM, the case 7 = 4/3 could be a nice toy model of a matter dominated universe approaching near dust conditions, 
though it would be a better model than that provided by 7 = 5/3 as it would include the MWCB. 



C. The drawbacks. 



It emerges from the discussion in the previous subsections that the k = 1 models, whether 7 = 4/3 or 7 = 5/3, 
share two unappealing features: (A) The maximal and minimal values of the rest mass density {mc^n or mc^n^) 
respectively occur along the symmetry centers r = and r = tt, while the maximal and minimal values of the internal 
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energy density (p — mc^n « 3p/2 or p — m(P"n} « 3p) occur along r = tt and r = 0. This means that along r = 
all energy density is rest mass energy density and there is a higher density of particles but zero temperature and 
pressure for all times (see discussion further ahead), while at ?■ = tt there are less particles by volume but a higher 
temperature and pressure. (B) The fact that pressure and temperature vanish along r = for all evolution times 
(forced by the specific spacial dependence of the metric contained in the function F) is a very different situation from 
p, T tending to zero asymptoticaly as the fluid expands and cools. This situation leads to the entropy per particle 
(18c) diverging as r 0, so that the left hand side of the Gibbs equation (14) (Ts'), remains finite (as the right hand 
side does). Although this feature is problematic, it is consistent with the fact that an absolute zero temperature (if 
T is stricly zero instead of T ^ asymptoticaly) is impossible to attain in equilibrium classical thermodynamics [ p4[ . 
This feature (appearing also for the cases A: = 0, — 1) is perhaps the strongest drawback on the physical plausibility of 
the models. 



D. Final remarks. 

The models derived and presented in this paper are mathematicaly simple, theoreticaly consistent and satisfy 
physicaly motivated equations of state. Whether one considers it a virtue or a defect, the models bear an extremely 
close resemblance (observationaly and dynamicaly) to FLRW cosmologies in the matter dominated regime with a 
dust source, therefore they share the possible failure of matter dominated standard cosmology in dealing with recent 
observational evidence pst . However, even if hard core cosmologists might dismiss them as not sufficiently realistic, 
the solutions can at least be applied as self-consistent cosmological toy models or as useful educational tools in teaching 
cosmology. Hopefuly the results of the present paper, besides rescueing the Stephani Universes from oblivion, will 
motivate and improve further research based on this nice and elegant class of solutions. 
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FIG. 1. Age of universe today. This figure displays Hoito — ti) vs fio for the = 1 model with 7 = 5/3 (solid line) and 
for a fc = 1 FLRW cosmology with a dust source. The values of Ho{to — ti) follow from (58) with A = 1 and Ai = 10~^, with 
Ho ~ 60 km (mpc seg)~^. Notice the almost identical fit between both cosmological models. The same close fit (not displayed) 
with standard cosmology occurs for k = —1, as well as for the k — 1, k — —1 models with 7 — 4/3. For the case A; = standard 
cosmolgy yields Qq — 1, while for the Stephani models (either 7 = 4/3, 5/3) we would obtain a plot analogous to the case fc = 1 
with f2o > 1- These (and the remaining) plots were obtained with the help of the symbolic computing package Maple V. 



FIG. 2. The relative ratio of the differences of Ho{to — ti) for the Stephani and FLRW models. This plot displays 
the ratio Xq = logj^g 1^) ~ ^01/^0 as a function of ilo for the case fc = 1 of 7 = 4/3, 5/3, where Yo = Hoito — ti) and Yq is the 
FLRW equivalent. The smallness of this ratio illustrates the close fit in the dynamics of the Stephani models and standard 
cosmology. The resemblance is tighter for 7 = 4/3 (curve below) than for the model 7 = 5/3 (curve above). 
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